Growth of spatial correlations in the aging of a simple structural glass 
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We present a detailed numerical study of dynamical heterogeneities in the aging regime of a 
simple binary Lennard- Jones glass former. For most waiting times and final times t, both the 
dynamical susceptibility X4(ii iw) and the dynamical correlation length ^4{t, t^) can be approximated 
as products of two factors: i) a waiting time dependent scale that grows as a power of t-u,, and ii) a 
scaling function dependent on t, t^ only through the value of the intermediate scattering function 
C{t,tw)- We find that Xiitjtw) is determined only in part by the correlation volume. 



PACS numbers: 64.70.Q-, 61.20.Lc, 61.43.Fs 

When a liquid's temperature is rapidly reduced, it can 
become supercooled and eventually undergo a transition 
into a glass state. As the transition is approached, the 
relaxation time and the viscosity grow by several orders 
of magnitude. The glass transition is the point at which 
the relaxation time of the liquid becomes longer than a 
fixed laboratory timescale^. Consequently, a system in 
the glass state is out of thermodynamic equilibrium. In 
particular, physical aging is observed: time translation 
invariance (TTI) is broken, i.e. correlations between ob- 
servablcs at times (the waiting time) and t (the final 
time) depend non-trivially on both times, and not just 
on their difference t — t^"^. 

Several observations about the glass transition, includ- 
ing the rapidly increasing relaxation timescales, the pres- 
ence of non-exponential relaxation, and the violation of 
Stokes-Einstein relations between viscosity and diffusiv- 
ity, have motivated the assumption that dynamical het- 
erogeneities are present in the supercooled liquid"'^''^'^. 
Dynamical heterogeneities are nanometer-scale regions 
containing molecules that rearrange cooperatively at very 
different rates compared to the bulk'^'^. More recently, 
dynamical heterogeneities have been directly observed 
in local probe experiments in supercooled liquids^'^, 
glasses and granular systems near the jamming tran- 
sition^''. However, except for the pioneering work by 
Parisi^^, most simulations studying dynamical hetero- 
geneities in off-lattice models of structural glasses have 
focused on the supercooled regime^^'^^'^*'^^'^^'^^'^^'^^ 

In this Letter we present a simulation study of dynamic 
spatial correlations of a continuous-space, quasi-realistic 
glass model in the aging regime, with the goal of char- 
acterizing their extent, space dependence, and time de- 
pendence. We discuss the spatial correlations of density 
fiuctuations, which have been used before to characterize 
dynamical heterogeneities in the supercooled regime^"^'^^. 
We first present detailed results from an isotropic probe, 
which allows us to focus on the strength and spatial ex- 
tent of the correlations, and later on briefly discuss the 
case of an anisotropic probe, which permits studying the 
geometry of the correlated region. A theoretical frame- 
work that postulates the presence of local fiuctuations 
in the age of the sample^^'^^, predicts that probabihty 
distributions of local fiuctuations in the aging regime are 



approximately independent of t^ at fixed values of the 
two-time global correlation Cgiobai(i, t^,). This prediction 
has been confirmed in simulations of both spin giasses^^ 
and structural glasses^^ (where Cgiobai is the self part 
of the intermediate scattering function). Motivated by 
this scaling, we examine the possible presence of scaling 
behavior as a function of Cgiobai and t^, in the dynamic 
spatial correlations. 

In order to probe the spatial density correlations, we 
use the 4-point (2-time, 2-position) correlation fimc- 



tion 



13. 



g4{r,t,t^) = -^(^"^Sir ~ rkitu.) + r,{t^)) (1) 

^ ik 

xwi\r,{t) - r,(tj| V(|r,(t) - r^it^)])) - (^^)'; 

where Q{t,t^) = Ylf^i w{\ri{t^)-ri{t)\) and io(|ri-r2|) 
is an overlap function which is unity for |ri — r2| < a 
and zero otherwise, (a = 0.3 cr^^ is an upper bound 
for the typical amplitude of the vibrational motion^"^.) 
By integrating g4{r,t,tw) over space, the generalized 4- 
point density susceptibility X4 = / d^^ 34 (r, t, t^)^^ is ob- 
tained. The self part of the intermediate scattering func- 
tion C = Cgiobai(i,i™) = ^ LiLi exp(iq.(rj(t)-rj (<„,))), 
can be used as a measure of the correlation between the 
configurations of the system at times t.^, and t: it is unity 
when t = tw and decays to zero when t — i^, iiu^'^'^°- 

We performed 4000 independent molecular dynamics 
runs for the binary Lennard- Jones (LJ) system of Ref.^°, 
which has a mode coupling critical temperature Tc = 
0.435. A system of 8000 particles was equilibrated at a 
temperature Tq = 5.0, then instantly quenched to T = 
0.4, and then allowed to evolve at that temperature for 
2.5 X 10'^ Lennard- Jones time units. The origin of times 
was taken at the instant of the quench. After the quench, 
data were taken at times in a geometric sequence t„ — 
tor"^, (with n = 0, . . . , 61, io = 10, r = lO^/^s « i_09). 

We first discuss the 4-point density susceptibility Xi- 
In the top panel of Fig. 1 we see that for fixed wait- 
ing time t.m, as the time difference t — t.^ increases, the 
spatial correlation increases till it reaches a maximum, 
and then decreases. The value of the correlation xa at 
the peak and the time difference t — t^ at which the 
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FIG. 1: (Color online) Top panel: Generalized density sus- 
ceptibility X4 as a function of the time difference t — t^, for 
waiting times tw ~ 13.8, • ■ ■ , 1318.26. Middle panel: Rescaled 
density susceptibility versus 1 — C, for the same waiting times. 
Bottom panel: Rescaling factor of the density susceptibility, 
X4,{C — 1/e), versus the waiting time. The straight line 
is a fit to X4(C' ~ 1/e) — at^^, with fitting parameters 
b = 0.353 ± 0.003 and a = 0.83 ± 0.02. The curved line is 
a logarithmic fit: ^4,{C = 1/e) — c\n{tjn) + d, with parameters 
c = 1.57 ± 0.04 and d = -2.6 ± 0.2. Statistical errors are 
shown by error bars in the data for all three panels. 



peak occurs both grow as a function of t^. Analogous 
behaviors are observed in numerical simulations of su- 
percooled liquids (see, e.g., Ref.^'^'^^), as the temper- 
ature is decreased; and in experiments in 2D granular 
systems^°, as the area fraction is increased. As we an- 
ticipated above, we test for a possible scaling behav- 
ior with C = Cgiobai by plotting, in the middle panel 
of Fig. 1, the rescaled quantity X4/X4(C = 1/e) as 
a function of 1 — C We observe that all curves ap- 
proximately collapse onto a single master curve. Thus 
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FIG. 2: (Color online) Top panel inset: Rescaled S4,{q,t,tw) 
as a function of q for two (t,„ , t) pairs, with low-q fits shown 
(dashed lines). The correlation lengths extracted from the 
fits are indicated in the key. Top panel: Correlation length 
as a function of t — t^ for various tw's. The "+" represent 
the {t — tw,^4) values for the two curves in the inset. Middle 
panel: Rescaled ^4 versus 1 — C. Bottom panel: Rescaling 
factor of the correlation length, £,4{C = 1/e), versus tw The 
straight line is a fit to £,4,{C = 1/e) = tt(i„)'', with parameters 
a = 0.853 ± 0.007 and b = 0.126 ± 0.001. The other line is a 
logarithmic fit: ^4(C = 1/e) — c\n{tw) + d, with parameters 
c = 0.195 ± 0.001 and d = 0.651 ± 0.008. Statistical errors are 
shown by error bars in the data for all three panels. 

we can say that Xii^^tw) ~ xt{tw)4'iC{t,tiu)), where 
Xl(tw) = Xiit,tw)\c{t,t^)=i/ei a rescaling factor that de- 
pends only on tw, and 4>{C{t,tw)) is a scaling function 
which depends on times only through the value of the 
intermediate scattering function C{t,ty^). As we see in 
the bottom panel of Fig. 1, the rescaling factor as a 
function of t^, can be fitted with the power law form 
Xl{t^) = atj, with the power b = 0.353 ± 0.003^5. As 
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FIG. 3: (Color online) 4-point density susceptibility rescaling 
factor versus the dynamic correlation length rescaling factor. 
The straight line corresponds to 0(^4)'', with b = 2.89 ± 0.03 
and a = 1.23 ± 0.02. Statistical errors for both abscissas and 
ordinates are shown for all the data points. 



the figure shows, this fit is much better than a fit to a 
logaritlimic form cln(tu,) + d. 

We now extract a dynamic correlation length ^4{t,tw) 
from our data. By Fourier transforming the correla- 
tion function gi(r,t,tw) we obtain the 4-point dynamic 
structure factor S'4(q, i, t^,)^^'^^'^^. We fit its small q 
behavior with an empirical scaling form^^: S'4(q) = 
{S2 - K) fiUH) + K, with fix) = 1/(1 + x''). Choos- 
ing 7 = 2 would give an Ornstein-Zernicke (OZ) form, 
but there is theoretical evidence^"* that OZ may not be 
a good description of S'4(q), and in fact 7 = 2.9 ± 0.1 
provides substantially better fits to our data^''. Since 
the q-dependent part of 5'4(q) scales with £,4\q\, and its 
value is reduced from its maximum at |q| = to half of its 
maximum at |q| = 1/^4, we interpret ^4 as a correlation 
length. The fits are performed in the range 0.4 < q < 1.9, 
for 175 curves corresponding to all the time pairs shown 
in Fig. 2. The fitting parameters S^it^t^), £,4{t,tuj) and 
K{t, tw) are determined independently for each time pair 
(i, t^). The scaling function f{x) is not allowed to change 
with {t,tm): the (t, tu,)-independent empirical parameter 
7 is determined by minimizing the average square error 
over the whole set of fits. The maximum difference be- 
tween the data and the corresponding fits never exceeds 
7.5%^^. Two typical fits are shown in the inset of the top 
panel of Fig. 2. 

The top panel of Fig. 2 shows the correlation length 
^4 as a function of the time difference t — for = 
13.8, 1318.26. We see two regimes in this plot: i) an 
early aging regime^^ {t — tw ^ t^), where ^4 is approxi- 
mately time translation invariant (i.e. it depends on t~t^ 
but not on t^), and ii) a late aging regime [t — t^ }^ tyj), 
where ^4 generally grows with increasing t^ at fixed t—ty^. 
This increase of ^4 with t^ is analogous to its increase in 
supercooled Hquids as the temperature is lowered^^'^*. 

We find that .^4 always grows as a function oi t — t.^, 
at fixed waiting time tw Simulations in the supercooled 
regime show divergent results on this issue: in Ref.^^ 



the correlation length increases initially, reaches a maxi- 
mum and goes to zero for long time differences, while in 
Refs.^^'^® the correlation length always grows with time 
difference. Similar to X4, the correlation length ^4 also 
displays an approximate scaling behavior with Cgiobai: 
in the middle panel of Fig. 2, we plot ^4/^4(C ~ 1/e) 
as a function of 1 — C. The curves collapse onto a sin- 
gle master curve and we can conclude that ^4(t,tu,) ~ 
£.4{tw)'p{C{t,t^)), where CK^^-) = U{t,tw)\c{t,t^,)=i/e- 
In the bottom panel of Fig. 2, we observe that the de- 
pendence of the scaling factor ^Kt^) on can be fit 

roughly equally well by a power law form ^| {tw) = a{tw)'' 
where b = 0.126 ± 0.001^^; or by a logarithmic form 
CK^iu) = cln(tu,) + d. However the power law form seems 
preferable because it is consistent with the physically mo- 
tivated scaling xl ~ {^4)^ found below. 

To test whether the value of Xiit^tw) is controlled by 
the correlation volume (^4(t, i^,))'^, in Fig. 3 we plot the 
susceptibility scaling factor xlitw) as a function of the 
correlation length scaling factor ^|(tm), and find that 
the data are indeed well represented by a power law form 
with a power close to 3 : X4° « ($4°)'' with b = 2.89±0.03 
(see also Ref.^^). However, looking at the data in more 
detail reveals an additional variation of Xiitjtw) that 
cannot be due to the correlation volume: if we compare 
the two plots of Xiit^tw) and £,4{t,tyj) versus the global 
correlation C{t,t^) at constant t^ (the middle panels of 
Figs. 1 and 2), we find that as C{t,tu,) — > 0, Xiit^'tw) 
almost reaches zero, but i,4{t,tyj) increases. Clearly, the 
evolution of X4 at constant t^, is controlled not only by 
the correlation length but also by the amplitude of the 
correlations. 

We now briefly turn our attention to the anisotropic 
spatial structure of the correlations^^'^*'^^. We de- 
fine a new spatial correlator ^/4(r, k, i, tu,) by replac- 
ing, in the definition of g4{v,t,tyj), all occurrences of 
w[\r[t) — r(t„)|) with the anisotropic function cos(k • 
(r(i) - r(t^)))23. We choose k = (6.68,0.00,2.67), with 
|k| = 7.2, i.e. at the peak of the static structure func- 
tion S{k). The Fourier transform of t/4(r, k, t^,) with 
respect to r is the 4-point anisotropic dynamic struc- 
ture factor iS4(q, k, f^,)^^. As we see in the left panel 
of Fig. 4, this 4-point function is indeed anisotropic: 
the contours of constant 54 extend further in the qy 
direction than in the q^ direction. We fitted 1S4 with 
an empirical form which generalizes the isotropic one: 
54 = [Si -K)f ([(g,6,,)2 + [qyU.vfY'^) + with 
f{x) = 1/(1 + x'^). Here 52,i^,^4,x, and ^4,;^ are {t,t^)- 
dependent fitting parameters, and 7 = 3.4 ± 0.1 (in- 
dependent of (t, i^,)) is determined by a global fit to 
all the data shown in Fig. 4. The middle panel of 
Fig. 4 shows that the rescaled quantity (^(q, k, i, tu,) = 
(54(q, k,t,i^) - K{t,t^))/{S2{t,t^) - K{t,t^)), com- 
puted for six time pairs, collapses onto one curve that 
agrees well with f{x). The right panel of Fig. 4 shows 
a contour plot of (T(q, k, !„,). Here the contours for 
various time pairs collapse, as long as a is kept constant 
and rescaled wavevector components S,4.xqx and S,4.yqy are 
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(Color online) Left panel: Contours for constant 54(q, k, t, t^,), plotted as a function of the and qy coordinates, for 
waiting times — 17.38, 100, 575.44, and constant Cgiobal — 0.3. The correlations extend more in the qy direction 



than in the qx direction. Middle panel: Rescaled 4-point function cr(q, k, t, t^) for 6 time pairs, for = 17.38, 100, 575.44, and 
Cgiobal = 0-3) 0-5- -A^U of the rescaled data approximately collapse onto one curve which is fairly well fitted by the scaling form 
1/(1 + x''). Right panel: Contours of constant CT(q, k, t, tm) = 0.2,0.45,0.7, plotted as a function of the rescaled wavevector 
i^A:xqx,^4,yqy), for three time pairs; = 17.38, 100, 575.44 with Cgiobal = 0-3. 



used on the axes. Since ky = 0, and q- = 0, 54(q, k, t, t^u) 
is probing correlations between displacements in the x di- 
rection. We find that ^4.a;(i,itu) > S,4:.y{t,tw), i.e. the 
correlation extends further along the direction of the 
probed displacements than along the perpendicular di- 
rection, and we interpret this as a manifestation of co- 
operative "string-like" motion, like in the supercooled 
regimc^^'^^'^^. 

In summary, we have studied the aging regime of a 
simple glass-forming model and found evidence of scaling 
for the dynamical susceptibility Xii^j '^w) and the dynam- 
ical correlation length ^4(t,tu,) as products of a power 
of tw times a scaling function that depends on times 
only through Cgiobal (t, iu>). Additionally, we have found 
evidence of time-independent scaling functions that de- 



scribe the q dependence of both the isotropic and the 
anisotropic dynamic structure factors. In the anisotropic 
case, we have found evidence for "string-like motion" , 
similar to the one observed in supercooled liquids. 
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